Quasiclassical theory of twin boundaries in High-T c superconductors 
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We investigate the electronic structure of twin bound- 
aries in orthorhombically distorted high-T c materials using 
the quasiclassical theory of superconductivity. At low tem- 
peratures we find a local instability to a time-reversal sym- 
metry breaking state at the twin boundary. This state yields 
spontaneous currents along the twin boundary that are mi- 
croscopically explained by the structure of the quasiparticle 
bound states. We calculate the local density of states and find 
a splitting in the zero-bias peak once the spontaneous currents 
set in. This splitting is measurable by STM techniques and is 
a unique signature of time-reversal breaking states at a twin 
boundary. 

A large number of experiments have convincingly 
demonstrated that Cooper pairs have basically d x 2_ 
wave symmetry in high-temperature superconductors 
In a strict sense this classification applies only to super- 
conductors with perfectly tetragonal crystal symmetry 
where the d x 2_ y 2 represents a "relative angular momen- 
tum" of the Cooper pairs with reduced symmetry. On 
the other hand, there is a number of slightly orthorhom- 
bically distorted systems such as YBa2Cu307 (YBCO) 
where the intrinsic crystal deformation removes the sym- 
metrical distinction between this pairing state and a con- 
ventional s-wave type. We may interpret this also in the 
way that the orthorhombic distortion couples the two 
pairing channels (|]|]. Clear evidence for this kind of 
mixing has been found in recent c-axis Josephson exper- 
iments between YBCO and Pb These experiments 
established the presence of an s-wave component and the 
dominance of the d x 2_ y 2-wave symmetry. 

An interesting aspect of the s-wave admixture due to 
orthorhombic distortion occurs in the vicinity of twin 
boundaries (TB) which separate the two degenerate or- 
thorhombic lattice shapes (twins). It was suggested that 
the two order parameter components, s- and d-wave, 
could locally twist in a way that time reversal-symmetry 
T is broken. Various physical properties are connected 
with this effect. In Ref. Q it was shown on the level of 
Ginzburg-Landau theory that the T- violating state is ac- 
companied by a spontaneous current flowing parallel to 
the twin boundary. Recently this system was studied by 
solving the Bogolubov-de Gennes (BdG) equations m 
confirming the results of the Ginzburg-Landau descrip- 
tion and discussing the local /^-characteristics observ- 
able in a tunneling experiment. Time- reversal breaking 
surface states were also predicted [Q and experimentally 
found @ at surfaces of YBCO. 

In this paper, we study this type of system and its 
properties using the quasiclassical theory of supercon- 



ductivity developed by Eilenberger ||. We calculate the 
structure of the order parameter self-consistently at all 
temperatures. Our calculation includes the orthorhombic 
distortion in form of anisotropic quasiparticle masses. In 
contrast to Ref. the twin boundary has no extension 
and is not explicitly pair breaking. In the bulk this for- 
mulation gives a mixing between s- and d-wave pairing 
due to the orthorhombic distortion in agreement with the 
other methods. For the twin boundary we show that the 
time-reversal symmetry is indeed broken at low enough 
temperature. The lower symmetry is noticeable also in 
the energy levels of quasiparticle bound states at the TB. 
In the T-invariant state there is a large density of bound 
states at zero-energy. The breakdown of time-reversal 
symmetry in the TB leads to the splitting of this energy 
level, and this effect in the local density of states is ob- 
servable by STM (scanning tunneling microscopy). The 
properties of these bound states lead also to a microscopic 
interpretation of the spontaneous current along the TB. 

In YBCO the orthorhombic symmetry reduction is par- 
tially induced by lattice deformation and mainly due to 
the presence of CuO-chains. The introduction of this or- 
thorhombicity into the quasiclassical formulation is done 
by changing the single electron dispersion to 
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where the (dimensionless) constant — 1 < c < 1 parame- 
terizes the distortion (the sign of c defines the two twin 
domains). Here the x- and y-axes correspond to the di- 
rection [1,1,0] and [1,-1,0] of the tetragonal lattice. This 
leads to a non-cylindrical FS with hp (6) = v / 1 fc e F s ° in2 fl 
where 9 is the angle between k and the A^-direction and 
/cfo the Fermi wave vector of the undistorted system. We 
will neglect all effects of electron motion in c-axis direc- 
tion. With this anisotropy the averaging over the FS has 
to be modified compared to the standard case, 
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Here, Nq is the normal density of states at the Fermi 
energy, and 
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where m is a normalization factor. A realistic estimate 
for the parameter c is ~ 0.2; this is mainly due to the 
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chains [Q. Besides this modification of the dispersion 
we will neglect all other effects related to the chains. 

In the following, we restrict our discussion to the s- 
wave and d x 2_ y 2-wave pairing channels. The pairing in- 
teraction is modeled as 



V{9, 9') = V s + 2V d sin(26>) sin(20') , 
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where the angular dependence of the d-wave part follows 
from our choice of the coordinate frame. The coupling 
strengths V s /a are eliminated in favor of the bare critical 
temperatures T cs /d and the energy cut-off cj c , the latter 
assumed to be the same for both channels. Then the 
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coupling constants can be defined as V d 



n>0 (" + V2) 

The bulk solution for the off-diagonal quasiclassical 
Green's function is f u {9) = A(0)/f2(0), where Q(9) = 
y^u> 2 + \A(9)\ 2 . The gap function is a linear combination 
of the s- and d- wave component, A(0) = A^ sin(20) +A S . 

The self-consistent calculation leads to two coupled 
equations for the s- and d-wave component. In the case 
c = we either find a one-component solution (i.e., the 
s- or the d-wave component vanishes), or both are finite 
and appear in the complex combination d ± is, i.e. they 
have the relative phase = ±7r/2 Q. At T = the 
complex phase occurs for 0.6 < T cs /T c d < 1.0 (see the 
inset in Fig. |l|). 

An expansion in c immediately shows that the s and 
d components are coupled for c ^ and a finite value of 
one of the components drives the other component to be 
non-zero. A non- vanishing c also leads to a renormalized 
onset temperature for superconductivity, 
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where c = c£^ 27rTc l/(n + 1/2). 

The numerical results for the two order parameters and 
their relative phase at T — are shown in Fig. |l| as a 
function of T cs /T c d- They always coexist. In contrast to 
the case c = 0, the relative phase cj> varies continuously 
and is different from zero in a narrow window of values 
of T cs /T c d- The region of broken time- reversal symmetry 
is shrinking with increasing c|, in qualitative agreement 
with the GL approach. 
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FIG. 1. Zero-temperature properties of an orthorhombic 
superconductor with mixed d- and s-pairing. For different 
magnitudes of the orthorhombic distortion c = 0.2 (main 
graph) and c = (inset), the two order parameters and the 
relative phases are shown as a function of T C3 /T c d- 



In the following we will assume T cs /T c d = 0.3 (see 
Ref. 0) and c = 0.2 which leads to A s ~ 0.25A d and <j> = 
or 7r (see Fig. [|), i.e., no bulk T- violating state occurs 
in agreement with the present experimental status. 

We will now consider a TB, a boundary between 
two domains with the relative orientation of the crys- 
tal axes of 90°. The quasiclassical equations have to 
be solved along classical trajectories, which are char- 
acterized by a momentum direction on the Fermi sur- 
face fcp At the TB the k- vectors with the same k v - 
component have to be matched. The reflection coefficient 
for a given fcp in the absence of a boundary potential is 
R = («Fxi. _ «Fx2) 2 /(wfx1 + w Fx2) 2 - If we consider the 
energy dispersion given in Eq. (nf) (c = sign(a;)|c|), we 
find that wpxi — ^Fx2 for all scattering directions. Con- 
sequently we have no normal reflection at the TB B . 




FIG. 2. Fermi surfaces and bulk gap (shaded region) on 
both sides of the TB. Fermi wave vectors fcf-vectors and 
Fermi velocities vf for a perpendicular trajectory are indi- 
cated by arrows. The parallel components of vf for these 
trajectories on both sides of the TB have different sign. 



We solve the quasiclassical equations numerically in a 
self-consistent way for various temperatures for a single 
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twin boundary. The calculation was simplified by ap- 
plying the Schopohl-Maki transformation to the Eilen- 
berger equations p2[ . Consistent with experimental ob- 
servations we assumed that the d-wave component 
is identical on both sides, whereas the s-wave component 
changes sign. 
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FIG. 3. Order parameters and relative phase in the vicin- 
ity of a twin boundary for T cs /T c d = 0.3 and c = 0.2. Upper 
graph: T- invariant state of the TB for T = 0.1T cd . The 
s-component goes to zero at the TB. Lower graph: locally 
T- violating state for T = 0.01T ct j. The absolute value of 
the s-component is finite for x = and the relative phase 
changes smoothly from to n. Inset upper graph: Modulus 
squared of (imaginary) s-component as a function of tempera- 
ture. Inset lower graph: spontaneous current density in units 
of W^eNoT^vp. 



ticular feature of local T-violation is the appearance of 
spontaneous currents in the vicinity of the TB flowing 
in opposite directions on both sides. These spontaneous 
currents at T = 0.01T C£ j are shown in Fig. || (inset lower 
graph). Note that screening effects are not taken into 
account here. 

The local density of states (LDOS) can be calculated 
by solving the Eilenberger equation for real energies. The 
results are given in Fig. ^. In the T-invariant state a 
pronounced zero-energy peak occurs in the vicinity of 
the TB. In the T-violating state, the zero-energy peak 
splits into two symmetric peaks. The separation of the 
peaks is proportional to A s (0). The LDOS can be mea- 
sured by scanning tunneling microscopy (STM). In par- 
ticular, the splitting of the zero-energy peak would prove 
the existence of a time-reversal breaking state in YBCO 
||. The temperature at which the splitting occurs, 
T ~ 0.05T C£ ; ~ 0(5K), is well within reach of current 
low-temperature STM technology. 



N(E,x) 




E/T cd 



N(E,x) 




"iE/T cd 



Results for two different temperatures are shown in 
Fig. [| For the higher temperature T = 0.1T cd - the rela- 
tive phase jumps from to 7r at the twin boundary and 
A s goes to zero directly at the TB. On the other hand, 
for T = 0.0lT c d we observe a local T-violation where the 
relative phase changes continuously from to ir with a 
non- vanishing s-wave pairing component. The inset in 
the upper graph shows |A S | 2 at the TB as a function of 
temperature. It becomes finite for T < 0.05. A par- 



FIG. 4. Local density of states in the vicinity of a twin 
boundary for T C3 /T c d = 0.3 and c = 0.2. The two graphs 
correspond to T = 0.1T ct j (upper graph) and T = 0.01T c< j 
(lower graph). In the T-invariant state there is a large peak at 
zero energy. The peak splits symmetrically in the T-violating 
state. 



To investigate the microscopic nature of the gap fea- 
ture in the LDOS and the spontaneous current it is useful 
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to turn to an angle-resolved view of the problem. The 
zero-energy peak in the LDOS in the T-invariant situ- 
ation can be understood by considering classical trajec- 
tories with angles 9 centered around 0, tt and ±7r/2 in 
a range 289 = arcsin(A s /A<;). The corresponding An- 
dreev (electron-hole) bound state for such angles is con- 
fined between two superconducting domains with a phase 
shift of 7r which leads to a zero-energy level jl3| . On the 
other hand, broken time-reversal symmetry will gener- 
ate phase shifts different from n leading to up or down 
shifted energy levels and the splitting of the zero-energy 
peak. Alternatively this effect may be also seen as a 
driving mechanism: the large DOS of the zero-energy 
bound states at the Fermi energy gives rise to a local 
Fermi surface instability opening a (pseudo) gap. 

How does this splitting lead to spontaneous currents? 
Each state (bound or extended) carries a current parallel 
to VF(k). The phase gradient seen by these states leads 
to a shift in energy (bound states) or spectral weight (ex- 
tended states). As a result left- and right-moving states 
are differently occupied. At T = only the right-moving 
of the bound states are occupied. On the other hand, 
due to the shift of the spectral weight the left-moving 
continuum states give the dominant contribution of the 
continuum states. They cancel the perpendicular current 
contribution of the bound states (as required by current 
conservation), leaving only a parallel component of the 
total current Q] . The sign of the parallel current is dif- 
ferent on both sides of the TB due to the symmetry of 
the dispersion relation (the momenta and velocities for 
one of these states are indicated in Fig. ||). 

In conclusion, we have studied the electronic structure 
of twin boundaries in orthorhombically distorted high- 
T c superconductors. The orthorhombic distortion was 
introduced in the quasiclassical formalism of supercon- 
ductivity and the corresponding equations were solved 
self-consistently. In contrast to Ref. || we do not as- 
sume that the TB has a pair-breaking effect. At low 
temperatures, we found a localized T-violating state at 
the TB. These localized states create spontaneous quasi- 
particle currents that flow in parallel to the TB. In addi- 
tion, we find a splitting of the zero-bias anomaly in the 
local density of states. Using realistic parameters we es- 
timate that this splitting could be observed around 5K, 
i.e., using available low-temperature scanning tunneling 
microscopes. This splitting would be a unique signature 
of T-violation in such systems. 
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